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Starting from the experimental evidence that high-energy nucleus-nucleus collisions cannot be 

described in terms of superpositions of elementary nucleon-nucleon interactions, we analyze the 

possibility that memory effects and long-range forces imply a nonextensive statistical regime during 

high energy heavy ion collisions. The relevance of these statistical effects and their compatibility 

with the available experimental data are discussed. In particular we show that theoretical estimates, 

obtained in the framework of the generalized nonextensive thermostatistics, can reproduce the shape 

r^ ' of the pion transverse mass spectrum and explain the different physical origin of the transverse 

O^ \ momentum correlation function of the pions emitted during the central Pb+Pb and during the p+p 

^\ ■ collisions at 158 A GeV. 
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I. INTRODUCTION 



In the last years, many efforts have been focussed on the study of high energy nuclear collisions, resulting in a 
better understanding of the strongly interacting matter at high energy density. Large acceptance detectors provide 
^ ' a detailed analysis of the multiplicity of produced particles; this in turn can be related to possible signatures of the 
f^ formation of a new phase of matter, the so-called quark-gluon plasma (QGP), in the early stages of high energy 
t^ heavy ion collisions ||l|. 

One of the most interesting experimental results is that high energy nucleus-nucleus (A+A) collisions cannot be 

described in terms of superpositions of elementary nucleon-nucleon (N+N) interactions (proton-proton or proton- 

^^ \ antiproton). The different conditions of energy density and temperature at the early stage of the Pb+Pb and of 

^^ the p+p collisions (to mention just two extreme examples) generate collective effects that modify the features of 

freeze-out observables. Of course the most appealing explanation would be to interpret the presence of these evident 

experimental differences as an indirect consequence of the formation of QGP in the early stages of heavy ion collisions. 

It is an experimentally established fact H that the slope parameter of the transverse momentum distribution 

y ■ significantly increases with the particle mass when going from p+p to central Pb+Pb collisions. On the other hand, 

^ ' the preliminary results of the NA49 Collaboration [plH] indicate that the transverse momentum fluctuations in central 

r* Pb+Pb collisions at 158 A GeV have a quite different physical origin with respect to the corresponding fluctuations 

. 1^ , in p+p collisions. Let us also recall that the observed J/^ suppression increases continuously and monotonically 

^^ ■ from the lightest (p+p) to the heaviest (S+U) interacting nuclei, but it exhibits a clear departure from this " normal" 

5h ' behavior for central Pb+Pb collisions [p| . 

These experimental features have been explained in terms of collective effects in the hadronic medium (such as the 
presence of a transverse hydrodynamical expansion) and considering the strong influence of secondary rescatterings 
in heavy ion collisions [g||2-|§]. It is a very important issue to understand whether the system is able to reach full 
thermalization; up to now this question has been the subject of many experimental and theoretical studies. However 
several experimental observables, such as transverse mass spectrum, multiplicity of particles, three-dimensional phase- 
space density of pions, are well reproduced in the framework of local thermal equilibrium in a hydrodynamical 
expanding environment [jSHl^. 

In this paper we start from the hypothesis that memory effects and long-range forces can occur during high energy 
heavy ion collisions: in this case the strong influence of collective effects on the experimental observables can be 
understood in a very natural way in the framework of the generalized nonextensive thermostatistics. In this context 
we show that the transverse momentum correlations of the pions emitted from central Pb+Pb collisions can be very 
well reproduced by means of very small deviations from the standard equilibrium extensive statistics. 



'Accepted for publication in Eur. Phys. J. C 



In SecO we summarize the main assumptions contained in the derivation of the dynamical kinetical equations and 
their relevance to the determination of the equilibrium phase-space distribution. In Sec. Ill we examine whether these 



conditions are met in the high energy heavy ion collisions and which experimental observables can be sensitive to the 



presence of nonextensive statistical effects. In Sec. IV we introduce the nonextensive Tsallis thermostatistics which is 
considered the natural generalization of the standard classical and quantum statistics when memory effects and/or 
long range forces are not negligible. In Sec.[V| we show how the generalized distribution function modifies the shape of 



the transverse mass spectrum and in Sec.Vl we investigate the influence of nonextensive particle fluctuations in the 



determination of the measure of the transverse momentum correlations. Finally, a discussion of the results and some 



conclusions are reported in Sec. VII, 



II. DYNAMICAL KINETICAL APPROACH TO THERMAL EQUILIBRIUM AND EXTENSIVE 

STATISTICS 

Assuming local thermal equilibrium, the freeze-out observables are usually calculated within the extensive thermo- 
statistics. In particular, the equilibrium transverse momentum distribution is assumed to be, in the classical case, the 
standard Maxwell-Boltzmann (Jiittern) distribution. If the emitted particles are light enough, such as pions, to have 
noticeable quantum degeneracy, one must use the Bose-Einstein (or Fermi-Dirac, for fermions) distribution. 

When the system approaches equilibrium, the phase-space distribution should be derived as a stationary state of 
the dynamical kinetical evolution equation. If, for sake of simplicity, we limit our discussion to the classical case, the 
Maxwellian distribution is obtained as a steady state solution of the Boltzmann equation. Let us now briefly review 
the principal assumptions in the derivation of the Boltzmann equation and examine how these approximations can 
affect the determination of the equilibrium distribution p3| . One important assumption concerns the collision time, 
which must be much smaller than the mean time between collisions. This request can be expressed by the condition 
nrQ <IC 1, where n is the density and rg is the effective range of the interactions. This condition has two important 
physical consequences, a) There is no overlapping between subsequent collisions involving a given particle and the 
interactions can be described as a succession of simple binary collisions, b) It is always possible to define a time 
interval in which the single particle distribution does not change appreciably and its rate of change at time t depends 
only on its instantaneous value and not on its previous history. Hence this property reflects the Markovian character 
of the Boltzmann equation: no memory is taken into account. 

The second assumption is analoguous to the first one, but for the space dependence of the distribution function: 
its rate of change at a spatial point depends only on the neighborhood of that point. In other words the range of the 
interactions is short with respect to the characteristic spatial dimension of the system. 

The last important assumption is the so-called Bolzmann's Stosszahlansatz: the momenta of two particles at the 
same spatial point are not correlated and the corresponding two body correlation function can be factorized as a 
product of two single particle distributions. This assumption is very important because it allows us to write down 
the collisional integral (which is equal to the total time rate of change of the distribution function) in terms of the 
single particle distribution only. The above assumptions, namely the absence of non-Markovian memory effects, the 
absence of long-range interactions and negligible local correlations, together with the Boltzmann's H theorem (based 
on the extensive definition of the entropy) lead us to the well-known stationary Maxwellian distribution. 

The basic assumption of standard statistical mechanics is that the system under consideration can be subdivided 
into a set of non-overlapping subsystems. As a consequence the Boltzmann-Gibbs entropy is extensive in the sense 
that the total entropy of two independent subsystems is the sum of their entropies. If memory effects and long range 
forces are present, this property is no longer valid and the entropy, which is a measure of the information about the 
particle distribution in the states available to the system, is not an extensive quantity. 

In the next section, we will see whether the above assumptions are implemented during high energy nuclear collisions 
and if there appear experimental signals that can be interpreted as a consequence of the presence of a nonextensive 
regime. 

HI. DO HEAVY ION COLLISIONS SATISFY EXTENSIVE STATISTICS? 

It is a rather common opinion that, because of the extreme conditions of density and temperature in ultrarelativistic 
heavy ion collisions, memory effects and long-range color interactions give rise to the presence of non-Markovian 
processes in the kinetic equation affecting the thermalization process toward equilibrium as well as the standard 
equilibrium distribution p4Ml7|. 



A rigorous determination of the conditions that produce a nonextensive behavior, due to memory effects and/or 
long-range interactions, should be based on microscopic calculations relative to the parton plasma originated during 
the high energy collisions. At this stage we limit ourselves to consider the problem from a qualitative point of view 
on the basis of the existing theoretical calculations and experimental evidences. 

Under the hypothesis that QGP is generated in the high energy collisions, the quantities characterizing the plasma, 
such as lifetime and damping rates of quasiparticles, are usually calculated within finite temperature perturbative 
QCD. This approach is warranted by the fact that at high temperature (beyond a few hundred MeV) the strong 
QCD coupling ag = g^/47r becomes very small and "weak coupling" regime takes place. In this regime {g <^ 1) the 
longitudinal gluon propagator is characterized by a Debye mass m^ « gT, and the corresponding screening length 
(Aj-, = m~^ w l/.9^) is much larger than the mean interparticle distance ((r) !v n'-^'^ « ^/T) (for a review, see 
e.g., Ref. [p8[). Therefore a great parton number is contained in the Debye sphere and the ordinary mean field 
approximation of QGP (Debye-Hiickel theory) holds. Nevertheless, in the proximity of the phase transition, the 
QGP is no longer a system of weakly interacting particles and non-perturbative QCD calculations become important. 
Recently it has been shown |19| that for temperature near the critical one (T > l.lTc) non-perturbative calculations 
imply an effective quark mass very close to T, sensibly different from the value, proportional to gT, obtained in 
perturbative regime. Near the phase transition the characteristic quantity nvQ « "-Af) should then be very close to 
one and only a small number of partons is present in the Debye sphere: the ordinary mean field approximation of 
the plasma is no longer correct and memory effects are not negligible. In addition, we observe that in high density 
quark matter the color magnetic field remains unscreened (in leading order) and long-range color magnetic interaction 
should be present at all temperatures. 

From the above considerations it appears reasonable that, if the dcconfining phase transition takes place, non- 
Markovian effects and long range interactions can influence the dynamical evolution of the generated fireball toward 
the freeze-out stage. Moreover they will affect the equilibrium phase-space distribution function if thermal equilibrium 
is attained. A signature of these effects should show up in physical observables. In this context, we notice that the 
authors of Ref. M raise a controversy on the Markovian description of the multiple scattering processes, one of the 
assumptions of the UrQMD model. This example concerns ( e~^e~ —^ fj,'^ fi~) scattering, where a determination of the 
angular distribution based on the Markovian approximation of the transport theory leads to wrong results pfl,pT^,0| . 

The aim of the present work is to suggest a possible interpretation of the above mentioned difference between 
the correlactions/fluctuations of the pion transverse momentum measured in p+p and Pb+Pb collisions [O as a 
signature of the nonextensivity of the system. In fact, it has been shown p2] that, whether the equilibrium condition 
is realized or not, the measure of the correlations/fluctuations strongly depends on the factorization of the multiparticle 
distributions and on the event by event multiplicity of the particles under consideration. Such property implies that 
the (multi)particle distribution functions are not inclusive (the distribution functions of one, two, or more bodies, 
depend on the presence of the other particles of the system) . In turn, this is a manifestation of a nonextensive behavior 
of the system. 

In addition we notice that, from a recent analysis of the average pion phase-space density at freeze-out in S- 
nucleus and Pb-Pb collisions at SPS, the experimental data indicate a slower decrease with increasing p± than the 
Bose-Einstein curve [ p3[ . The analysis of the same quantity in 7r-p collisions seems, instead, to be consistent with 
the standard expectations. We remind that in Ref. [g^ the authors show that by assuming the standard Jiittern 
momentum distribution of emitted particles in high energy collisions one is led to unrealistic consequences and the 
physical freeze-out is actually not realized. 

From the foregoing considerations we conclude that both theoretical calculations and experimental observables 
agree with the existence of nonextensive features in high energy heavy ion collisions. 

IV. GENERALIZED NONEXTENSIVE TSALLIS STATISTICS 

Several new developments in statistical mechanics have shown that in the presence of long-range forces and/or 
in irreversible processes related to microscopic long-time memory effects, the extensive thermodynamics, based on 
the conventional Boltzmann-Gibbs thermostatistics, is no longer correct and, consequently, the equilibrium particle 
distribution functions can show different shapes from the conventional well known distributions. Standard sums, or 
integrals, which appear in the calculation of thermostatistical quantities, like the partition function, the entropy, the 
internal energy, etc. can diverge. These difhculties are well known in several physical domains |E5[ since long time. 

A quite interesting generalization of the conventional Boltzmann-Gibbs statistics has been recently proposed by 
Tsallis p9| and proves to be able to overcome the shortcomings of the conventional statistical mechanics in many 



physical problems, where the presence of long-range interactions, long-range microscopic memory, or fractal space- 
time constraints hinders the usual statistical assumptions. Among other applications, we quote astrophysical self- 
gravitating systems [g^ , the solar neutrino problem [g^ , cosmology [g9| , many-body, dynamical linear response theory 
and variational methods [Q, phase shift analyses for the pion-nucleus scattering | |3l[ |. 

The Tsallis generalized thermostatistics is based upon the following generalization of the entropy [E6| 
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where pk is the probability of a given microstate among W different ones and g is a fixed real parameter. 

The new entropy has the usual properties of positivity, equiprobability, concavity and irreversibility, preserves 
the whole mathematical structure of thermodynamics (Legendre transformations) and reduces to the conventional 
Boltzmann-Gibbs entropy S = — '^^Pi logpi in the limit q ^ 1. 

The deformation parameter q measures the degree of nonextensivity of the theory. In fact, if we have two indepen- 
dent systems A e B, such that the probability oi A + B is factorized into pa+b{ua, ub) = Pa(ua)Pb{ub), the global 
entropy is not simply the sum of their entropies but it is easy to verify that 

Sg{A + B) = Sg{A) + Sg{B) + (1 - q)Sg{A)Sg{B) . (2) 

For a better appreciation of the meaning of nonextensivity in a physical system, we discuss here another important 
property. Let us suppose that the set of W microstates is arbitrarily separated into two subsets having Wl and Wm 
microstates {Wl + Wm = W), respectively, and define as pi, = ^iJlPi and pm = Si^vF +iP« ^^"^ corresponding 
probabilities: hence pl +Pm — ^- It is then easy to show that 

Sq{{Pl}) = Sq{pL,PM) +pI Sq{{p^/pL}) + p\j Sq{{p^/pM}) , (3) 

where the sets {pi/pL} and {pi/pn} are the conditional probabilities. This is a generalization of the famous Shannon's 
property but for the appearance of p|^ and p^^ instead of p^ and pM in the second and third terms of the right hand 
side of (Q). Since the probabilities {pi} are normalized, pf > pi for q < 1 and pf < pi for q > 1: as a consequence 
values oiq<l{q>l) will favour rare ( frequent) events, respectively. 

The single particle distribution function is obtained through the usual procedure of maximizing the Tsallis entropy 
under the constraints of keeping constant the average internal energy and the average number of particles. For a 
dilute gas of particles and/or for q w 1 values, the average occupational number can be written in a simple analytical 
form 1^ 

- ' (4) 
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where the + sign is for fermions, the — for bosons and /? — 1/T. In the limit q —^ 1 (extensive statistics), one recovers 
the conventional Fermi-Dirac and Bose-Einstein distribution. 

Let us remind that the Tsallis generalized statistics does not entail a violation of the Pauli exclusion principle 
and does not modify the inclusive behavior of the bosons, but it modifies, with the generalized nonextensive entropy 
(|l|), the extensive nature of the standard statistics. At the equilibrium, the nonextensive statistics implies a finite 
temperature particle distribution different from the standard Fermi-Dirac and Bose-Einstein distributions; in the 
classical limit, one has the following generalized Maxwell-Boltzmann distribution [p6| : 

{n,), = [l + {q-l)f5{E,-^i)]^/^^-'^K (5) 

When the entropic q parameter is smaller than 1, the distributions (Q) and (^ have a natural high energy cut-off: 
Ei < l/[/3(l — q)] + A*, which implies that the energy tail is depleted; when q is greater than 1, the cut-off is absent and 
the energy tail of the particle distribution (for fermions and bosons) is enhanced. Hence the nonextensive statistics 
entails a sensible difference of the particle distribution shape in the high energy region with respect to the standard 
statistics. This property plays an important role in the interpretation of the physical observables, as it will be shown 
in the following. 



V. TRANSVERSE MASS SPECTRUM AND Q-BLUE SHIFT 

The applicability of the equilibrium statistical mechanics relies on the recognition that the experimental transverse 
mass spectrum and the relative abundance of charged particles (pions, kaons, protons, etc.) can be described within 
the equilibrium formalism. The single particle spectrum can be expressed as an integral over a freeze-out hypersurface 
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where g is the degeneracy factor and f{x,p) is the phase-space distribution. 

Indeed, by assuming a purely thermal source with a Boltzmann distribution, the transverse mass spectrum can be 
expressed as follows: 

HN 

^A m^Ki (z) , (7) 



m±dm 



where z = rn±/T and Ki is the first order modified Bessel function. In the asymptotic limit, m± ^ T {z ^ 1), the 
above expression gives rise to the exponential shape 

dN , , 

= B V^^e-^ , (8) 

which is usually employed to fit the experimental transverse mass spectra and provides an indication of a thermal 
energy distribution. 

Although, from a general perspective, the experimental distribution of particle momenta in the transverse direction 
is described by an exponentially decreasing behavior, at higher energies the experimental data deviate from the usual 
Boltzmann slope. This discrepancy can be cured by introducing the dynamical effect of a collective transverse flow 
which predicts the so-called blue shift factor, i.e., an increase of the slope parameter T at large m\^. 

The experimental slope parameter measures the particle energy, which contains both thermal (random) and col- 
lective (mainly due to rescattering) contributions. The thermal motion determines the freeze-out temperature Tfoi 
namely the temperature when particles cease to interact with each other. In the presence of a collective transverse 
flow [pj one can extract Tfo from the empirical relation 

T^Tfo + m{vi_f , (9) 

where (wx) is a fit parameter which can be identified with the average collective flow velocity; the latter is consistent 
with zero vci p -\- p reactions, but can be large in collisions between heavy nuclei (where rescattering is expected to be 
important). In Eq.(^), m is the mass of the detected particle (e.g. pions, kaons, protons). 

Let us consider a different point of view and argue that the deviation from the Boltzmann slope at high pj_ can 
be ascribed to the presence of nonextensive statistical effects in the steady state distribution of the particle gas; the 
latter do not exclude the physical effect of a collective flow but rather incorporates a description of it from a slightly 
different statistical mechanics analysis. 



If long tail time memory and long-range interactions are present, as already discussed in Sec. Ill, the Maxwell- 
Boltzmann distribution must be replaced by the generalized distribution (g) . We consider here only small deviations 
from standard statistics (5 — 1 « 0); then at first order in (g — 1) the transverse mass spectrum can be written as 

-^^ ^Cm^ Ik, [z) + ^^^z^ [3 K,(z) + K:,{z)]\ , (10) 

where K3 is the modified Bessel function of the third order. In the asymptotic limit, z ^ 1 (but always for q values 
close to 1, such that (q — l)z <C 1), we obtain the following generalization of Eq.(g): 

dN ^ ,__ f , g - 1 2 
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z + :^^—z'\ . (11) 



From the above equation it is easy to see that at first order in (g — 1) the generalized slope parameter becomes the 
quantity Tq, defined as 



Tg^T+iq-l)m^. (12) 

Hence nonextensive statistics predicts, in a purely thermal source, a generalized q-blue shift factor at high m±; 
moreover this shift factor is not constant but increases (if g > 1) with m± — y^m^ + p^, where m is the mass of the 
detected particle. We have in this case an effect very similar to the one described by the phenomenological Eq.(0). We 
remark that the proposed distribution of Eq.(0), based on the nonextensive thermostatistics, has the same high p± 
dependence of the parameterization reported in Ref. 139] to fit the transverse distribution for central Pb+Pb collisions. 

The value of the entropic parameter g is a measure of the nonextensivity of the system and, as a consequence, it 
should be fixed by the reaction under consideration. On the basis of the above considerations, at a fixed reaction 
energy, one should find that larger value of q are required for heavier colliding nuclei, because at higher energy density 
we have a larger probability for the phase transition to be achieved. 

A quantitative description of the particle production in nonextensive statistics goes beyond of the scope of this 
paper; it is fair to assume that Eq.(0), for q > I value, can be in good agreement with experimental results in heavy 
ion collisions. In this framework, indeed, one can understand the slower decrease (for heavier colliding nuclei) of the 
pion phase-space distribution with increasing p± . 

It is worth noticing that for small deviations from the standard statistics (Kg — 1)| < 0.05) the experimental 
midrapidity transverse mass distributions are well reproduced with the same slope parameter of Ref. , within the 
statistical errors. 

In the next Section we shall calculate the measure of the transverse momentum fluctuations in the framework of 
the equilibrium nonextensive statistics. We notice that, since small nonextensive statistical effects are consistent with 
the experimental transverse mass spectrum, we can use, in what follows, the same temperature extracted from the 
standard analysis of the experimental results. 

VI. TRANSVERSE MOMENTUM FLUCTUATIONS IN NONEXTENSIVE STATISTICS 

In Ref. [ p4| , Gazdzicki and Mrowczyhski have introduced a quantity <&, which measures the event-by-event fluctu- 
ations and does not explicitly depend on the particle multiplicity, but appears to be sensitive to the correlations. It 
is defined as 



*^ = VW"^' (13) 

where a; is a variable such as energy or transverse momentum, z^ — x — x is a single particle variable and Zx = 
X]i=i(^i — x) is the corresponding sum-variable, TV being the number of the particles in the event. 

A non-vanishing value of the measure $ implies that there is an effective correlation among particles (dynamical 
or statistical) which alters the momentum distribution. It is usually assumed that if the experimental value of $ 
turns out to be constant in going from N+N to A+A collisions, then nucleus-nucleus collisions can be described as an 
incoherent superposition of elementary N+N scattering processes and the correlations among particles are unchanged 
p3,B3. We wish here to point out that this might not be the case: indeed it will be shown that similar values of <& 
can be attributed to different statistical correlations, which in turn reflect distinct dynamical situations. 

Recently the correlation measure $p^ of the pion transverse momentum has been measured by NA49 Collaboration. 
The experimental results correspond to a value $^^p = 0.6 ± 1 MeV for Pb+Pb at 158 A GeV and to ^^^ = 5+1 
MeV in the preliminary measurements of the p+p collisions at the same energy [y,^]. However, in their most recent 
analysis of the Pb+Pb data, the NA49 Collaboration has estimated a contribution of A<I>pj^ = 5 + 1.5 MeV from the 
statistical two-particle correlation function alone and an "anti-correlation" contribution of A<I>pj^ — —4 + 0.5 MeV, 
stemming from the limitation in the two-track resolution of the NA49 apparatus [Q . Hence it appears that the actual 
physical values both for p+p and for Pb+Pb collisions are nearly equal. In any case, as it was observed in the same 
paper, the physical origin of transverse momentum fluctuation remarkably changes from Pb+Pb to p+p collisions 
and, although the pion transverse correlation measure has been studied within many different theoretical approaches, 
the obtained results are somewhat controversial and not well understood |qj35|-|3^. 

We will show in the following that the physical differences in the origin of the pion transverse momentum fluctuations 
in Pb+Pb collisions with respect to the p+p ones can be understood in the framework of the nonextensive statistics. 
For this purpose, keeping in mind that the whole mathematical structure of the thermodynamical relations is preserved 
in the nonextensive statistics, it is easy to show that the two terms in the right hand side of Eq.(|3|) can be expressed 
in the following simple form 
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P±-P±] {n)q, (14) 



and 



{Z'pJ 1 /• d^p 



-pjwf{p--~^-y^''^'^^' ^''^ 



(N) 
where the quantity pj^ is the average transverse momentum, 

_ 1 f d^p , , f d^p , ^ 

In the above equations we have indicated with (n) q the mean occupation number of Eq. (0) extended to the continuum 
and with (An^)^ = {n'^)q — (n)'^ the generahzed particle fluctuations, given by 



where E is the relativistic energy E — ^vn? -V p^ . For q—\, one recovers the wefl-known text-book expression for 
the fluctuations of fermions (— ) and bosons (+). Eq.dlS) shows that the correlation measure $ arises from the particle 
fluctuations; this result does not come as a surprise since it is well known that the particle fluctuations are related to 
the two-body density or to the isothermal compressibility by means of the so-called Ornstein-Zernike relation |T3] . 

In the standard extensive statistics (g — 1), one obtains that the correlation measure $ is always negative for 
fermions and positive for bosons, while it vanishes in the classical standard statistics [p5|. Indeed (Z^) is defined 
in terms of the fluctuations An^, which, for the ideal Fermi-Dirac (Bose-Einstein) gas, are suppressed (enhanced) 
with respect to the classical case. In the nonextensive Tsallis statistics, the fluctuations of an ideal gas of fermions 
(bosons), expressed by Eq.(p7|), are still suppressed (enhanced) by the factor lq= (n)^; this effect, however, is modulated 
by the factor [1 + (g — l)/3(i? — [i)\~^ in the right hand side of Eq.(|l7[). Therefore the fluctuations turn out to be 
increased for q <\ (we remind that due to the energy cut-off condition the above factor is always positive) and are 
decreased for g > 1. The measure $ does not have, both for bosons and fermions, a constant sign and it is positive 
or negative depending on the value of the entropic parameter q. In the classical nonextensive case we also have 
(An^)q 7^ (?^)g- The quantity (An^)q is greater than (n)^ for g < 1, and smaller for g > 1 and the correlation measure 
$ can be different from zero for q ^ \. This is essentially due to the intrinsic nature of the system in the presence 
of memory effects and/or long-range forces and to the dependence of the (multi)particle distribution functions upon 
the surrounding medium. 

We use Eqs.(14) and (|l^) to evaluate the correlation measure $p^ for the pion gas (ttitt = 140 MeV), which we can 



compare with the experimental results measured in the central Pb+Pb collision by the NA49 collaboration |^,Q . Let 
us notice that the values of $p^ , which depends on T, [i and the entropic parameter g, are very sensitive to small 
variations of q. This is mainly due to the modified expression (|7|) of the fluctuations in the nonextensive statistics. 

In Fig. ^ wc show <I>p^ as a function of the parameter g, at different (fixed) values of the temperature and chemical 
potential. The experimental data are reported up to a maximum value of the transverse momentum pj_ = 1.5 GeV: 
hence we have limited the integration up to the upper limit of 1.5 GeV. Concerning the temperature T, we use the 
slope parameter derived from the fit of the transverse momentum spectra P,p3l. With T = 170 MeV and /i = 60 MeV 
(according to the suggestion of Q), we obtain <I>p^ = 5 MeV with q = 1.038, while the usual extensive statistics, 
with g = 1, gives $p^ = 24.7 MeV (we remind that a non-vanishing value of /i implies that there is no chemical 
equilibrium) . 

By choosing /i = MeV one can obtain the same value <I>p^ = 5 MeV with a smaller q {q — 1.015). In this case, 
with g = 1, $p^ = 13.6 MeV. 

Thus, a very little deviation from the standard Bose-Einstein distribution is sufficient to obtain theoretical estimates 
in agreement with the experimental results. We notice also that the value employed, g ~ 1, justifies, a posteriori, the 
approximation of the analytical form of the distribution function (Q) and the validity of Eq. (|l^) for the transverse mass 
spectrum. The correlations are much more affected by deviations from the extensive statistics than the single particle 
observables, such as the transverse mass distribution. We can obtain a very good agreement with the experimental 
results of $p^ , still being consistent with the other single particle measurements. 

In this context, it is worth noticing that, even for small deviations from the standard statistics, an appreciable 
difference in the distribution function and in the fluctuations (with respect to the standard ones) appears especially 



in the region of high transverse momenta. This behavior is consistent with the interpretation of the nonextensive 
statistics effects, since the detected high energy pions should be essentially the primary ones, hence more sensitive to 
the early stage of the collisions. In order to show this fact, we have evaluated the partial contributions to the quantity 
^P_L) by using Eq.s ( |l4| ) and (1^), and by extending the integration over p± to partial intervals Apx = 0.5 GeV. 
The results are shown in Fig.E, at T = 170 MeV and /i = 60 MeV. In the standard statistics (dashed line), $p^ is 
always positive and vanishes in the p^-intervals above w 1 GcV. In the nonextensive statistics (solid line), instead, 
the fluctuation measure ^p^ becomes negative for p± larger than 0.5 GeV and becomes vanishingly small only in 
px-intervals above ^ 3 GeV. 

In view of these considerations, we suggest an analysis of the Pb+Pb measurements which allows to disentangle 
the contributions to ^p^ in different pj_'s intervals. A negative value of ^p^ at high transverse momentum intervals 
could be a clear evidence of the presence of a nonextensive regime. 

VII. CONCLUSIONS 

The nonextensive statistics appears suitable to evaluate physical observables recently measured in heavy ion collision 
experiments. The spectrum of the transverse momentum can be reproduced by means of the nonextensive distribution, 
which naturally takes into account typical collective effects in heavy colliding nuclei such as the increasing of the slope 
parameter with highp^ and high particle masses. The calculated correlation measure ^pj^ agrees with the experimental 
value by considering a small deviation from the standard statistics. Within the nonextensive approach we have also 
found negative partial contributions to ^pj^ at high p± (larger than 0.5 GeV). In this regime the $pj^ value is not 
affected by resonance decays: hence an experimental confirmation of our prediction would be an unambiguous signal 
of the validity of nonextensive statistics in relativistic heavy ion collisions. 

The quantity q is not just a free parameter: it depends on the physical conditions generated in the reaction and, in 
principle, should be related to microscopic quantities (such as the mean interparticle interaction length, the screening 
length and the frequency collision into the parton plasma). We have found that a value of q slightly larger than one 



characterizes the correlation measure 4>, both by suppressing the boson fluctuations in Eq.(15) and by enhancing the 
importance of the region of large transverse momenta in Eq.dlJ). In the diffusional approximation, a value q > 1 
implies the presence of a superdiffusion among the constituent particles (the mean square displacement is not linearly 
proportional to time, having, instead, a power law behavior (x^) ex t", with a > 1). Effects of anomalous diffusion 
are strongly related to the presence of non-Markovian memory interactions and colored noise force in the Langevin 



equation |41|. In this context, we note that in Ref. [Q, the QCD chiral phase transition is analyzed in analogy with 
the metal-insulator one and it is shown that an anomalous diffusion regime can take place. The meaning of anomalous 
diffusion can be found in the presence of multiparticle rescattering, which is very large in Pb+Pb collisions nM ■ 

The early stage of the p+p collision is believed to be radically different from the Pb+Pb case and the influence of 
the nonextensive statistics appears not negligible in the last reaction. In these terms the anomalous medium effects 
and the remarkable experimental differences between light and heavy ion collisions should be better understood. 

Finally we notice that, for a given value of the deformation parameter q, the greater is the mass of the detected 
particle, the more strongly modified arc the fluctuations. For this reason, the present calculation predicts more 
important nonextensive effects for the correlations of heavier mesons and baryons, than for light particles: future 
measurements of ^p^ for these heavier particles should give a more clear evidence of this effect. Furthermore, the 
nonextensive statistics could imply a sizable modification of the hadronic mass compressibility, since this quantity is 
strongly related to particle fluctuations (being more sizably suppressed for particles of heavier masses^]) [Q . We also 
remind that the nonextensive statistics implies, for g > 1, a suppression of the standard boson fluctuations given in 
Eq. ( [l7|) , which is more effective at high transverse momentum whereas it appears negligible for low momenta. Shuryak 
suggests pO| that enhanced fluctuations for low p± bins, in the pion p± histograms, are a possible manifestation of 
disoriented chiral condensation. We may argue that this effect could be much more evident if one assumes the validity 
of the generalized statistics. New data and investigations are necessary to gain a deeper understanding of the high 
energy heavy-ion observables. 



^Using the same value, q = 1.02, of the deformation parameter, we can estimate a reduction of the standard Bose-enhancement 
factor 1 + (n) of about 2% for pions, of 7% for kaons and of 11% for p mesons. 
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FIG. 1. The correlation measure $p^ [MeV] as a function of the entropic parameter q, calculated in four different conditions 
of temperature T and chemical potential fj, (T = 170 MeV and ^ = 60 MeV solid line, T = 140 MeV and fj. = 60 MeV dashed 
line, T = 170 MeV and fj. = MeV dot-dashed line, T = 140 MeV and ^ = MeV dotted line). The value $p^(q = 1) 
corresponds to the correlation measure calculated within the standard extensive statistics. 
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FIG. 2. The partial contributions to the correlation measure "l>p^ [MeV] in different p± intervals, at T = 170 MeV and 
/i = 60 MeV. The dashed line refers to standard statistical calculations with q = 1, the solid line corresponds to 5 = 1.038 (note 
that by summing the contributions up to pj_ — 1.5 GeV we obtain the value $p^ = 5 MeV, indicated by the NA49 experiment 
M. Negative contributions, at high px steps, are predicted only within the nonextensive statistics. 
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